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ABSTRACT 

We perform a study of the spatial and kinematical distribution of young open clusters 
in the solar neighborhood, discerning between bound clusters and transient stellar con- 
densations within our sample. Then, we discriminate between Gould Belt (GB) a nd 
local Galactic disk (LGD) members, using a previous estimate (|EHas et al.ll2006ah of 
the structural parameters of both systems obtained from a sample of 0-B6 Hipparcos 
stars. Single membership probabilities of the clusters are also calculated in the sepa- 
ration process. Using this classified sample we analyze the spatial structure and the 
kinematic behavior of the cluster system in the GB. The two star formation regions 
that dominate and give the GB its characteristic inclined shape show a striking dif- 
ference in their content of star clusters: while Ori OBI is richly populated by open 
clusters, not a single one can be found within the boundaries of Sco 0B2. This is 
mirrored in the velocity space, translating again into an abundance of clusters in the 
region of the kinematic space populated by the members of Ori OBI, and a marginal 
number of them associated to Sco 0B2. We interpret all these differences by character- 
izing the Orion region as a cluster complex typically surrounded by a stellar halo, and 
the Sco-Cen region as an OB association in the outskirts of the complex. In the light 
of these results we study the nature of the GB with respect to the optical segment of 
the Orion Arm, and we propose that the different content of star clusters, the different 
heights over the Galactic plane and the different residual velocities of Ori OBI and 
Sco 0B2 can be explained in terms of their relative position to the density maximum 
of the Local Arm in the solar neighborhood. Although morphologically intriguing, the 
Gould Belt appears to be the result of our local and biased view of a larger star cluster 
complex in the Local Arm, that could be explained by the internal dynamics of the 
Galactic disk. 

Key words: (Galaxy:) solar neighborhood — open clusters and associations: general 
— (stars:) early-type — formation. 



1 INTRODUCTION 



The G ould Belt (GB) was firs t disco vered by John lHerschell 
l|l847h and Benjamin iGouldl (|l879l ') as a system of bright 
stars inclined with respect to the plane of the Milky Way. 
For more than a century, many studies have been devoted to 
describing its structure and its kinematical behavior, as well 
as to proposing a reliable global scenario that would account 
for its origin (for an extensive review on the subject, see 
Poppel 1997, 2001 and Grenier 2004). Today it is considered 
that, in the sco pe of the most recent theories of hierarchical 
star formation (|Efremovlll978l , Il995l : lElmegreen et al.ll2000l : 
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lElmegreenll2006l ). the GB is likely to be our closest giant 
star forming complex (jComeron 11200 ih . 

The stellar component of this complex takes the shape 
of a planar distribution of bright and youn g OB stars 
inclined with respect to the Galact ic plane l|Leshl Il968l : 
IStothers fc Frogell ll974l : IWestinI [l985l ) . Most of the young 
OB associations in the solar neighborhood are known to be 
part of the GB (Blaauw 1965; de Zeeuw et al. 1999; Elias et 
al. 2006a, hereafter Paper I). Also, a system of young, low- 
mass stars, detected by cross-ma tching X-ray and optical 
Hipparcos jPerrvman et al.| Il997^ based catalogs, appears 
to be associated with the GB ( Guillout et al.lll998l ). 



As we should expect from a giant ongoing star forming 
complex, the local interstellar medium is prominently asso- 
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ciated with the GB. The works by van den Berghl (|l966l ) on 
reflection nebulae , by Sandgvistl (11977) on dark clouds and, 
more recently, bv lGaustad fc Van BurenI |l993') on maps of 
infrared emission, have found a spatial distribution of the 
dark clouds of interstellar dust compatible w ith the pattern 
show n by the stellar component of the GB. IXomital ll 19861. 
Il987h and the study of the CO molecule by Taylor et al.l 
l|l987l ) seem to confirm this inclination for the local molec- 
ular clouds. 

Neutral hydrogen in the solar neighborhood has also 
been rela ted to the GB af t er the discovery of Lin dblad's "fea- 
ture A" l|Lindbladlll967l : iLindblad et al.lll973l ) . interpreted 
as a ring of gas with an expansion movement l|Qlandll982l : 
lElmegreenll 19821 ). After Lindblad's work, the giant molecu- 
lar clouds were foun d to be related to the most prominent 
OB star assoc iations (|Sancisi et al.|[l973 : lKutner et al.ll 19771 : 
Ide GeudlTggj ). A full map of the CO molecule over the sky 
later confirmed that most of the se clouds within 1 kpc from 
the Sun follow the GB pattern (|Dame et al.ll 19931 ). 

Then, if the GB is a star formation complex composed of 
both young stars and associations and interstellar material, 
we would expect to find a population of young open clusters 
following the GB pattern. This is obvious if we consider the 
concept of young star cluster in its widest sense, i.e., repre- 
senting the dense in ner regions of the hierarchical structure 
of young star fields l|Elmegreenll2006l ). Nonetheless, in this 
study we want to distinguish between a young cluster under- 
stood as a mere stellar condensation, and a gravitationally 
bound system (that is, as a condensation that "has had suf- 
ficient time and gravitational self-attraction to get mixed 
by stellar orbital motions" , as put by Elmegreen 2006) . Our 
work focuses on the analysis of the cluster system, and on 
how these objects distribute and behave in comparison with 
the GB defined by the massive stars. 

Since the first systematic studies tha t led to a discrim- 
ination between th e GB and the LCD l|Stothers fc Frogel 
Il974l : IWestinlll985l ) there has been a great improvement in 
the number of cataloged open clusters, as well as in the avail- 
ability of their astrometric and physical data. Surprisingly, it 
is not possible to find in the scientific literature any work de- 
voted to the study of the open clusters membership to the 
GB and their distribution within this system before 2006. 
Thus, for many years, it has been accepted that the GB 
did not contain a s ignificant population o f bound clusters. 
Only very recently IPiskunov et al.l l|2006l ). in their analy- 
sis of the Galactic open cluster population, discovered an 
open cluster complex (OCC) that they associate with the 
GB. Although they find this OCC as a density peak within 
the spatial distribution of clusters in the solar neighborhood, 
they estimate OCC membership probabilities by kinematical 
methods, through the analysis of the tangential velocities. 

In two previous papers (Paper I; Elias et al. 2006b, here- 
after Paper II) we have studied the spatial distribution and 
the kinematic properties of the OB stars and associations in 
the GB. Our line of work thus leads in a logical way to the 
study of the distribution of young open clusters in the GB 
and their kinematic properties. Our analysis will be centered 
in the comparison of the GB morphology as obtained from 
the distribution of massive stars and clusters. This will rep- 
resent another step in the understanding of the nature of the 
GB, and will also contribute to the knowledge of how star 



formation mechanisms proceed to the formation of stellar 
clusters. 



2 OPEN CLUSTER SYSTEM 

2.1 Associations, stellar condensations and bound 
clusters 

Prior to describing the selection criteria for our observational 
sample, we want to punctuate some ideas about the concepts 
of hierarchical star formation, association, loose groupings 
and clusters as gravitationally bound physical systems. 

Stars are born from molecular gas clouds whose inter- 
nal structure can be characterized by a fractal dimension 
value that apparently ranges from 2.5 to 2.7 (e.g. Sanchez 
et al. 2005, 2007). Simulations of the collapse of gas clouds 
hint at a hierarchical structure of the ste llar formation, wit h 
clusters present in several dense cores (|Walsh et al.ll2006l ). 
During the past two decades, a set of observations has been 
collected which indicates that young stellar groupings show 
hierarchical patterns that range from the larger scales of fioc- 
culent spiral arms and star complexes to the smaller scales 
of OB associations, OB subgroup s, small lo ose groups, clus- 
ters and cluster subclumps ( Efrcr novll 1995h . In other words, 
the newly born stars seem to follow the same fractal pattern 
as the gas clouds from which they were originated. 

The largest scale of stellar grouping, the great star com- 
plexes, would be associated with the gas superclouds with 
masses the order of 10^ solar masses. According to some au- 
thors (e.g., Comeron 2001, Efremov & Elmegreen 1998), the 
Gould Belt, with its spatial scale the order of 1 kpc, would 
be the star complex closest to the Sun, and thus would come 
from a single gas cloud with a mass of a few million solar 
masses. 

What happens at smaller scales? We already know that 
most of the OB associations in the solar neighborhood, with 
typical sizes of about 80 pc, are mainly distributed along the 
plane of the GB. They represent the observable spatial scale 
immediately below the star complex. But, does it make any 
sense to talk about a typical scale within a fractal structure? 
The answer should be sought in the observational bias that 
is introduced by the age limit of the sample: OB associations 
are detected and selected as concentrations of OB stars with 
a typical age of about 10 Myr. The existence of some general 
correlation between t he duration of the star format ion and 
the size of the region (|Larsonlll98ll : lEimegreenll2006l ) implies 
that, for a typical age of 10 Myr, the typical size of the r egion 
is 80 pc l|Efremov fc Elmegreenlll998l : lElmegreenll2006l ). 

Star clusters are formed in the cores of giant molecular 
clouds; they represent the stellar groupings associated to the 
inner and densest regions of the gas, and can be interpreted 
as the unavoida ble result of star formation in hierarchically 
structured gas (|Elmegreenl |200(? ). However, only a few of 
these condensations will still be gravitationally bound after 
10 Myr. It has been estimated that 90% of the clusters lose 
a high fraction of their stars in the first 10 Myr of their lives 
(Fall ct al. 200d). Thus, after 10 Myr of life, it is possible 
to find star forming regions that maintain a large number 
of star clusters, while ot hers like NGC 604 in M33 do not 
contain a single cluster (|Mafz- Apellani j 1200 ll ) . This seems 
to depend on variations of the mean density of the clouds; 
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those where the average density is low will form stellar con- 
centrations in the cores of the clouds, but they will not have 
enough binding energy to keep a bound cluster when the gas 
leaves. Thus, the clusters observed within a star-forming re- 
gion at a certain moment could be representatives of two 
distinct physical systems: either star clusters, gravitation- 
ally bound and able to survive galactic tidal forces, or a 
mere stellar condensation with a mean lifetime the order of 
10 million years or less. 

Thus the analysis of the distribution of star clusters 
(both bound systems and transient condensations) could 
give us information about the history of star formation in 
the GB, as well as about the physical conditions of the gas 
from which they were born. 

2.2 Cluster sample 

We extract our sample from the Catalogue of Open 
Cluster Data (COCD) a nd its Extension 1, compiled by 
iKharchenko et al.l (|2005al lbh. This catalogue has the advan- 
tage of homogeneity over other existing compilations, and 
since we are also interested in working with space velocities, 
the fact that the COCD had catalogued proper motion val- 
ues in the Hipparcos system and newly determined radial 
velocity data was decisive in our choice. 

As our aim is to study the cluster distribution in the 
GB, we establish distance and age limits in our sample. The 
GB system should be well contained within a heliocentric 
radius of 1 kpc (Stothers & Frogel 1974; Westin 1985; Paper 
I); a nd since its age ha s been estimated between 20 and 90 
Myr iTorra et al.ll2000l '). we only keep clusters younger than 
100 Myr. 

Thus, our first selection (that will be reduced after out- 
lier elimination, as we explain in the following section) is 
composed of 93 open clusters, 83 of which have complete 
kinematical information. We calculate for every cluster its 
Cartesian Galactic coordinates {X, Y, Z), where X is posi- 
tive in the direction of the Galactic center, Y in the direction 
of Galactic rotation and Z perpendicular to the Galactic 
plane so that they form a right-handed, orthogonal frame. 
We also calculate their respective space velocities, {U, V, 
W), for those clusters with radial velocity data. 

2.3 Gould Belt and Local Galactic Disk clusters 

In Paper I we developed a three-dimensional classification 
method that allowed us to separate the GB stars from the 
local Galactic disk (LCD) stars by purely spatial criteria. 
This method considered that the LCD and the GB could 
be described as a distribution along two intersecting planes; 
working with a sample of Htpparcos 0-B6 stars we obtained 
an estimation of the parameters that defined those planes. 
In Paper II, we used those parameters to classify strictly 
by spatial criteria a new sample of OB stars with kinematic 
data. Similarly, we now use our classification method and 
the parameters found in Paper I to separate the GB open 
clusters from those of the LCD field. We want to stress that 
this separation is done only by means of their spatial posi- 
tion in the three-dimensional space. 

As described in Paper I, this classification method leads 
also to the identification of outliers; i.e., those objects that 



are too far away from the mean planes and thus lie in re- 
gions of low density of probability in the sample space. We 
have found 9 outliers in our cluster sample; after their elim- 
ination the remaining sample consists in 84 clusters, 76 of 
them with full space velocity data (listed in table[T]). Further 
information on the detection and meaning of outliers may 
be found in Section 2 of Paper I and references therein. 

Finally, the individual GB membership probability (P) 
for each cluster is listed in table [T] In total, 40 clusters have 
a probability P > 50%, and thus, following the Bayes mini- 
mum error rate decision rule, we consider them to be mem- 
bers of the GB. 

2.4 Bound and unbound clusters 

This sample may contain both bound and unbound clusters, 
as we have previously seen (Section 2.1). A simple criterion 
to select those condensations with a high probability of being 
bound clusters is an age cut-off, if we consider that all the 
objects older than 10 Myr can be called bound clusters, in 
the sense that they have survived to the high rate of infant 
mortality that happens during the first 10 million years in 
the life of a cluster. Once this critical threshold has been 
surpassed the mean life of a cluster, although very dependent 
on the environmental conditions, is usually larger than 1 
Gyr. 

It is evident, though, that such a selection is just an 
approximation to the problem, and that to determine if a 
stellar system represents a bound cluster, we must compare 
its density with the tidal density at its position in the galaxy. 
However, we lack the complete information to perform this 
detailed analysis, and thus we must resort to empirical classi- 
fication criteria. The COCD catalog provides three variables 
for each cluster that may give us some additional informa- 
tion about what we understand as young stellar grouping 
and cluster; these variables are: cluster radius, core radius 
and age. 

In Figure 1, we represent the cluster sample in the clus- 
ter radius-age space. It is interesting to note that the clus- 
ters older than 10 Myr in our sample seem to group around 
an elongated strip with a positive slope and that only a 
few objects present a cluster radius that deviates from this 
mean distribution. The most notorious case is that of NGC 
2264, which has a radius close to 20 pc. This cluster seems 
to present a high degree of sub structure, spatial as well as 
kinematical (jFfiresz et al.l[200^ . that its central potential 
has not been able to erase. Thus we eliminate this object 
from the sample, for its properties appear to differ in some 
structural aspects from that of "classic" bound clusters. 

We have also drawn in Figure 1 an upper envelope of 
the main cluster distribution with ages older than 10 Myr 
and younger than 100 Myr. There seems to be a natural sep- 
aration between the main distribution of clusters older than 
10 Myr, and those that show a radius larger than expected 
for their age. Based on this apparent separation, we have 
traced this upper envelope of the main cluster distribution. 
The exact location of the line has been drawn by consen- 
sus of the authors. It simply aims to represent a dividing 
line between the largest concentration of "probable" bound 
clusters (represented as filled circles in the figure), and those 
that in spite of their age show a clearly distinct behavior. 
The objects located under this envelope line can be consid- 
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ered, from a conservative point of view, as highly probable 
bound clusters. 

Using the cluster radius and core radius data, we have 
estimated a pseudo-concentration parameter in the form 
log(Rci/Rco), that we represent in Figure 2 versus the cluster 
radius. Circled squares indicate those clusters that, accord- 
ing to the previously exposed criteria, can be considered as 
probable bound clusters. We note that most of these objects 
are distributed, as in Figure 1, along a straight line. How- 
ever, two of these objects (Platais 6 and NGS 2546) show a 

clear separation from the general tendency. 

Using the online tools at WEBDA database (iMermilliodI 

1 19951 ) we have recalculated the parameters of these two clus- 
ters from their photometric data. This analysis indicates 
that the photometric diagrams of Platais 6 show a good 
fitting for a distance modulus of 7 and log(t)~ 6.5. Simi- 
larly, a visual inspection of the color-magnitude diagram of 
NGC 2546, as cataloged by WEBDA, seems to indicate that 
it is a very young cluster (log(t) ^ 7) with a rich popula- 
tion of pre-main sequence stars. This diagram is similar to 
that shown by NGC 2362 and other young clusters located 
in the third Galactic quadrant (DolEado ct al. 2006, 2007), 
meaning that we are not facing an object almost 100 Myr 
old. Thus we consider these two objects as clusters with ages 
inferior to 10 Myr, and consequently they do not belong - 
according to our criteria- to the group of clusters with a 
high probability of being gravitationally bound. 

In table [1] we present the classification of the sample 
that follows from this reasoning. The last column shows an 
indicator of the probability of being a bound cluster, accord- 
ing to the criteria previously discussed (1 stands for those 
objects that are "probable" bound clusters, and 2 stands for 
those we consider as transient stellar condensations). 



3 ANALYSIS AND DISCUSSION 
3.1 Spatial distribution 

A two-dimensional projection of the spatial distribution of 
open clusters in the sky is shown in Figure 3 (top panel). 
There we see, as we commented in Section 1, that GB clus- 
ters (filled circles) are mostly concentrated towards southern 
Galactic latitudes, and that only three of them (P95 — IC 
4665, (423 = CoUinder 359, and (456 = Stephenson 1) clearly 
rise above the Galactic plane {b > 10 degrees). According 
to our criteria, [1395 and [1456 are probable bound clusters, 
while tt423 could rather be a transitory stellar condensation. 
We must note that these three objects are located in the 
first Galactic quadrant and that, apparently, are not related 
to any OB association. 

We also represent in this figure the OB stars used in our 
analysis of the spatial structure of the GB (Paper I). The 
associations Sco OB2, Ori OBI, Per OB2 and Lac OBI are 
classically thought to be components of the GB (e.g Olano 
1982), and studies that followed the Hipparcos mission sug- 
gest tha t Vel OB2, Tr 10 and CoUinder 121 also belong to 
the GB (|de Zeeuw et al.lll999l ). although their position close 
to the line of nodes where the GB intersects the Galactic 
disk adds quite some uncertainty to this membership assig- 
nation. In any case, it is evident that the Sco OB2 and Ori 
OBI associations respectively mark the South and North 



Galactic extremes of the inclined plane which best describes 
the stellar system known as the GB, and that in a certain 
sense the geometrical characterization of the GB is defined 
by the relative position of these two associations. 

Thus, we have marked with a red cross the stars belong- 
ing to th e Sco OB2 assoc i ation , following the coordinates as 
given bv lde Zeeuw et al] (|l999l ). that situate the complex in 
the range I = 290° - 360°, b = -10°- 30°and D = 100 - 
220 pc, and that we have tagged as Scorpius-Centaurus in 
the figure. In the same manner, we have selected the stars 
belonging to Ori OBI as those in the range I = 197° — 215° , 
6 = -12°- -26° and D = 300 - 550 pc (Ide Zeeuw et al.lll999l ). 
These stars appear in Figure 3 as green triangles, and the re- 
gion is tagged as Orion. Estimations of the age range of these 
associations indic ate that Sco OB2 m embers are between 5 
and 20 Myr old (ISartori et al.ll2003h . while the typical age 
of Ori O BI members is comprised between some 10 ^ yr and 
11 Myr (jBriceno et al.ll2005l : [Hernandez et al.ll2006l '). 

Hence we are facing two associa tions with sizes and age s 
relatively similar, and according to Ide Zeeuw et ahl |l99^, 
with a number of probable Sco 0B2 members larger than 
that of Ori OBI. With these data it sounds reasonable 
enough to think that the star cluster population related to 
Sco 0B2 should be larger than that related to Ori OBI. But 
what we observe in Figure 3A is exactly the opposite. There 
is not a single cluster within the frontiers of Sco OB2, while 
we have detected 11 objects associated to Ori OBI ((73, 
(174, tt75, t|76, tt77, t|80, P016, pOlS, P019, ^1020, P021), 
six of them being "probable" bound clusters (tt73=NGC 
1981,tt74=NGC 1976, tt75=NGC 1977, P018, ttl019, P020). 
This is better observed in Figure 3 (bottom panel), where 
we have represented only those clusters cataloged as "prob- 
able" bound ones, tagged with number 1 in the last column 
of table [T] This figure shows how almost the totality of the 
"probable" bound clusters associated to the GB are located 
in the Orion and the Puppis-Vela regions. Moreover, almost 
all of the other objects show no relationship with other OB 
associations in the GB. 

The same phenomenology can be observed in Figures 4 
and 5, where we have represented the GB cluster distribu- 
tion in the spatial projections XY and XZ. The LCD clus- 
ters tend to distribute uniformly across the Galactic plane in 
Figure 4, but that is not the case for GB open clusters. We 
observe how the GB's distribution is quite clumpy; many 
of the clusters are located in the region of Vela, and we 
specially note that the region of Orion (in the third quad- 
rant of the XY plane, and around the most negative values 
of Z in the Figure 5) presents an important concentration 
of members. Moreover, if we consider only clusters with a 
higher membership probability (75%, as noted in Figure 6), 
and thus eliminating most of the clusters around the line 
of nodes in which the GB plane intersects the LCD plane, 
the GB is practically reduced to the Orion region, as if that 
were the solely cluster population of the GB. 

This analysis shows as an evident fact that while the 
Ori OBI association is related with an important popula- 
tion of star clusters (be they transitory stellar condensations 
or gravitationally bound systems), the Sco 0B2 complex 
does not include a single star cluster within its frontiers, 
yet it is extrem ely rich in massive OB stars and pre-main 
sequence stars l|Sartori et al.|[2003h . In other words, inside 
an apparently single star formation complex as the GB, and 
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separated only by about 500 pc, we find two star forming 
regions, well defined both spatially and kinematically, that 
present two clearly distinct modes of star formation. Orion 
is an OB association that presents a large number of stellar 
condensations seen as clusters, some of which appear to be 
gravitationally bound. On the other hand, Sco-Cen seems to 
be forming only isolated stars or loose groups that do not 
present the shape of a star cluster, and thus have not been 
detected as such. But, how does this fact translate to the 
velocity space? 

3.2 Kinematic behavior 

It seems obvious that the different spatial position of the 
centroids of the two stellar groups analyzed in Section 3.1, 
implies also a different location in the velocity space, spe- 
cially in the V component due to their separation of almost 
500 pc in the X sods. The velocity centroids of these two 
associations, Ori OBI and Sco 0B2, as defined by the GB 
stars, are located at (-16.4, -9.5, -5.0) and (-8.0, -19.4, -6.0) 
km s~^, respectively. This can be clearly observed in Figure 
7 (top panel), where we have represented in the UV plane 
the isodensity lines of the GB as defined by massive stars, as 
well as the member st ars of Ori OBI a n d Sco OB2, accord- 
ing to the criteria by Ide Zeeuw et al.l (|l999l ) explained in 
Section 3.1, and the GB star clusters. Once more, as it was 
expected, we observe a clear correlation between the clus- 
ters' distribution and the Orion velocity centroid, while the 
number of GB clusters associated to the Sco-Cen complex 
is merely marginal. 

The difi'erence is even more evident if we limit ourselves 
to the youngest clusters, which could be representatives of 
transient stellar condensations (Figure 7, bottom panel). As 
we have commented above, the difference observed between 
the velocity centroids was predictable due to the differen- 
tial Galactic rotation and to the fact that both groups were 
quite separated in space. However, we wonder if the rotation 
field in this region of the Galaxy can completely explain the 
kinematic behavior of both the stars and the clusters in the 
GB.^ 

iMoreno et al] l| 19991 ) analyzed the velocity space of a 
sample of OB stars belonging to the GB, and found that 
the observed velocity field was not compatible with that 
obtained from star formation models, as well as with its 
dynamical evolution after the injection of momentum and 
energy in the primeval cloud. In other words, the residual 
veloc ities of th e stars showed a highly negative vertex devia- 
tion (|Filinlll95 7: Mihalas & Binncy 1981) that could not be 
explained by the dynamical model, because after a time in- 
terval quite inferior than the age of the GB, the differential 
rotation prevailed over the movements originated by the "ad 
hoc" energy and momentum inje ction, and gave place to a 
slightly positive vertex deviation. iMoreno et al.l |l99^ also 
found that if the stars belonging to the Sco-Cen association 
were eliminated, the vertex deviation became positive. That 
is, those models designed to explain the origin of the GB 
from supernovae explosions after a previous process of star 
formation, or from the impact of a high velocity cloud on the 
Galactic disk (see Poppel 1997, 2001 and Sartori et al. 2003 
for reviews about possibly formation mechanisms of the GB) 
are not able to explain the observed residual velocity field, 
unless the Sco-Cen stars are neglected. 



Now we want to evaluate the residual velocity ellipsoid 
for the GB clusters; in order to do so, we have corrected the 
velocities from solar motion and differential rotation using 
the Oort constants: A = 16 km s~^ kpc ~^ and B — -16 
km s"^ kpc ~^ (Paper II, Mendez et al. 2000). These cluster 
residual velocities are represented in Figure 8, along with 
the residual velocities of OB stars, where different marks in- 
dicate those belonging to Ori OBI and Sco OB2. Then we 
calculate the longitude of the vertex, for both systems, 
the LGD and the GB cluster samples. The result, if we cut at 
GB membership probability of 50%, is that Z„ = 9.3° ± 8.3° 
for the LGD, and that 1-, = -1.3° ±15.5° for the GB. This is 
a value very far from the GB vertex deviation of Z„ = —47° 
found in Paper II for the OB stars belonging to the GB. 
This is undoubtedly caused by the absence of clusters in 
the Sco-Cen association, that was responsible for the large 
vertex deviation of the GB (Moreno et al. 1999; Paper II). 
Moreover, if we keep only GB clusters with a membership 
probability higher than 75%, the vertex deviation of the sys- 
tem is Iv — 9.7° ±16.3°, which is practically the same as that 
of the LGD. 

The vertex estimation for the star cluster system in the 
GB thus gives us a double information. First, from a kine- 
matic point of view, the lack of star cluster associated to the 
Sco OB2 complex is demonstrated. Second, the difference in 
the velocity space between the Ori OBI and Sco 0B2 as- 
sociations cannot be completely explained by the Galactic 
differential rotation. 

As we may see in Figure 8, the residual velocity dis- 
tributions of these two associations show different behavior. 
While the stars belonging to Ori OBI present a main veloc- 
ity dispersion axis with a positive vertex, the Sco OB2 stars 
present a main axis that is almost perpendicular to the for- 
mer, with a clearly negative vertex deviation. Therefore, as 
it had already been noted (Moreno et al. 1999; Paper II), 
the vertex deviation in the solar neighborhood will depend 
on the selection of the sample. If most of the sample stars 
belong to Sco OB2, the vertex deviation will undoubtedly be 
negative. If we extend our sample farther away from the Sun 
in order to include Orion stars, we will reach some balance 
and thus the vertex deviation will turn towards values closer 
to zero. The latter is precisely what we observe in the star 
cluster population associated to the GB: the lack of clus- 
ters within the Sco-Cen complex makes the vertex deviation 
close to zero. 

3.3 Cluster complexes and scaled OB associations 

What we have observed when comparing the distribution of 
star clusters and OB stars in the GB is that the two great 
complexes that seem to define the North and South Galac- 
tic extremes of this large stellar structure, show a different 
behavior according to the scenery of hierarchical star forma- 
tion. Ori OBI shows a considerable portion of its stellar pop- 
ulation as grouped and forming star clusters, half of them 
being probable bound clusters with ages larger than 10 Myr. 
This kind of stellar system has been detected , observed and 
analyzed in both the Milky Way (Ef remov fc Sitnik lOSl; 
lAlfaro et al.lll99lh and other galaxies (e.g., 30 Doradus in 
the LMC; Hunter et al. 1995; Walborn et al. 2002); not only 
it contains a rich star cluster population, but also it is nor- 
mally immersed inside a stellar halo. 
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On the opposite side of the star forming regions concen- 
tration range, there can be found the scaled OB associations 
(SOBAs; e.g., Mai'z-Apellaniz 2001), of which NGC 604 in 
M33 is a good example. Although the star formation rate in 
NGC 604 is much higher than that observed in the Sco OB2 
association, they both have in common their lack of star 
clusters. Thus, separated only by 500 pc, there are two OB 
associations that, apparently, have been born from molecu- 
lar clouds under different ambient conditions, but that al- 
ways have been considered as the fundamental parts of a 
single stellar system known as the GB. However, Ori OBI 
seems to represent the stellar halo associated to a cluster 
complex, while Sco OB2 appears to be a clear example of 
an OB association, not related to cluster formation. 

If we interpret this result in terms of a hierarchical star 
formation process (|Elmegreenll2006l .[ 20081 ) ■ and considering 
that the age of the stars in both associations present similar 
ranges, we should consider either that the density maximum 
in the parental gas distribution of Sco 0B2 was inferior than 
the density peaks in the clouds that formed Ori OBI, or that 
the tidal forces in the Sco-Cen region were intense enough 
to destroy any substructure observable as a star cluster, in 
an interval of time lesser than 10 Myr. 

Any of these two possibilities requires some variations of 
the ambient physical conditions in spatial scales smaller than 
500 pc, be they due to an external difference of pressure that 
caused higher density peaks in the Orion region, or/and due 
to local gravitational potentials or shear forces that shorten 
the lifetimes of the transient stellar condensations in the 
Sco-Cen region. 



sent their residual velocities, corrected from solar motion 
and differential rotation. 

Ori OBI appears associated to the density maximum of 
the Orion Arm defined by the young clusters; on the con- 
trary, Sco 0B2 is located in the inner rim of the Arm, close 
to the Sun where the cluster density is lower. If we consider 
that a higher density of clusters is representative of a higher 
ambient pressure in the original gas, the relative position of 
these associations in respect to the main locus of the Orion 
spiral arm could explain their different content of star clus- 
ters. 

To better illustrate this, we define a Cluster Formation 
Index (CFI) that describes the relative content of clusters 
with respect to the OB stars that shape the GB: 

_ Cluster density 

Cluster density + OB star density 

The spatial densities of clusters and OB stars have been 
obtained by using gaussian kernels in the GB plane {X'Y'), 
with a a — 200 pc, and normalized so that the total sum of 
the density be equal to 1. Then, Figure 10 shows the distri- 
bution in the GB plane of the CFI parameter in the region 
that contains both OB associations. A line that joins the 
centroids of both associations has been represented over the 
density map and a cut along this line (Figure 11), clearly 
shows a CFI gradient that ranges from 0.62 at the maxi- 
mum near Orion to 0.39 in the vicinity of Scorpius Centau- 
rus. Since the distance to the Ori OBI centroid is a good 
estimator of the distance to the Orion spiral arm, this indi- 
cates that as we move away from the arm, the clustered star 
formation is less efficient. 



3.4 The Orion Arm 

Since the pioneering work of iBeckerl l| 19561 ). who traced the 
local spiral structure from the young star clusters distribu- 
tion within a radius of 2 kpc, each time that a new catalog 
of star clusters has been tailored, the corresponding map 
of this local spiral structure has been drawn again (e.g., 
Janes & Adler 1982). The results of doing this show a se- 
ries of cluster groupings with typical sizes the order of 1 
kpc l|Efremov fc Sitniki(l988l : lAlfaro et al.lll99ll , Il992l ) that 
seem to follow three segments of spiral arms which have re- 
ceived the names of Carina-Sagittarius Arm, Perseus Arm 
and Orion Arm (also known as Local Arm). The inclusion of 
these optical segments within the general scheme of the spi- 
ral structure of the Galaxy is controversial, and depending 
on the spiral tracers and the analysis techniques employed, 
different solutions for the number of arms, their pitch angle 
or the velocity of the density pattern have been found (see 
Naoz & Shaviv 2007, and references therein). 

In Figure 9 we represent the density map for the star 
clusters younger than 10 Myr, inside a square of side length 
4 kpc, centered in the Sun. The data have been extracted 
from the COCD catalog. In the figure we observe five main 
concentrations of young clusters that had previously been 
detected by other authors (e.g., Efremov & Sitnik 1988). In 
particular, the Orion and Cygnus complexes seem to align, 
delimiting the local optical segment of the Orion Arm. Su- 
perimposed over the young clusters density map, we have 
drawn the Ori OBI and Sco OB2 associations. The lines 
that depart from the centroids of both associations repre- 



3.5 The Nature of the Gould Belt 

The first detection of the GB was based in the fact that the 
brightest stars in the sky, specially in the southern hemi- 
sphere, present an eccentric position with respect to th e 
great circle of the Milky Way l|Herschelll 18471 : lGouldlll879l ). 
That is, it was a mere morphological matter. Later studies 
determined the main stellar and gaseous components of the 
GB, and from the analysis of their spatial and kinematical 
properties, the shape, size and kinematics of the GB were 
estimated. 

Although the main parameters that describe the GB, 
published along the past two decades, may differ depending 
on the authors (see Paper I; Poppel 1997 and references 
therein), the general scheme of the GB is similar for most of 
them: we are facing a star formation complex with a disk- 
like structure and a radius of about 500 pc, whose kinematic 
behavior is characterized by an expansion and a rotation 
with respect to an internal axis. In most of the works that 
contributed to the determination of these characteristics, the 
separation between the probable members of the GB and 
the stars belonging to the LCD has been performed. The 
analysis of the stellar component of both groups has led 
to the conclusion that the GB and the LCD are not only 
separated in the celestial sphere, but they present different 
kinematic properties as well. 

In most of these analyses it was necessary to perform 
a previous separation of the elements belonging to either 
group, and in all those cases such separation was mainly 
geometric (be it either in the l-b plane, in the different XY, 



Hierarchical Star Formation 7 



XZ, YZ projections, or via a three-dimensional analysis). In 
other words, the stars that are the brightest and the farthest 
from the Milky Way plane, and that belong to the Sco-Cen 
and Orion constellations, seem to have been the ones that 
opened the gates to this flood of studies, and also the ones 
that, still today, best define and delimit the geometry of the 
GB. 

Then, we shall perform the following experiment: let us 
consider the GB as defined only by two points whose coordi- 
nates in the phase space are given by the spatial coordinates 
of the centroids of the Ori OBI and Sco OB2 association, 
and their velocities by their respective central values of their 
residual velocities (table [2]). A schema of the geometry of the 
problem is shown in Figure 12. Using these values, we de- 
termine the inclination (i), the longitude of the ascending 
node {Q), the expansion velocity along the line that joins 
the spatial centroids (po) and the rotation velocity with re- 
spect to a point situated along the main axis of the system 
(loo), we reach the results listed on table [S] where we also 
present the range of values found in the literature for the 
different variables in the GB. 

These values are all very close to those of the fundamen- 
tal parameters of the GB that we can find in the literature 
(e.g.. Paper I; Fresneau et al. 1996; Lindblad 2000, among 
others). But on sight of Figure 9, where we represent the po- 
sition and the residual velocities of both associations with 
respect to the Orion Arm, we may ask ourselves if it makes 
any sense at all talking about a gas and stars system with 
spatial and kinematic coherence that may be described with 
some expansion and rotation velocities, or perhaps we should 
definitely drop the traditional hypothesis of a single, com- 
mon origin for all the features of the GB, and begin to look 
at it as a hazardous alignment -from out point of view- of at 
least two of the many clumps in the hierarchical structure of 
the Local Arm, with different densities and star formation 
history. In this sense, the GB would be simply the projec- 
tion over the sky of the recent star formation in the clouds 
close to the Sun but located far away from the fundamen- 
tal Galactic plane. Therefore, it is morphologically distinct, 
but does not necessarily represent a physical system with 
unique properties, different from the star formation within 
the Orion Arm. 

We believe that the distinct properties of these two 
clumps -their kinematic behavior, cluster content and height 
over the Galactic plane- may be explained by the internal 
dynamics of the Galactic disk. Further research should con- 
sider, as possible mechanisms involved in the development 
of such a structure, the passing of a density wave in a mag- 
netized medium or the presence of a long bar genera ting 
resonances in the external parts of the disk ((Gardner et al.l 
l2008l ). 

Thus, if we could observe our Galaxy from an external 
position several megaparsecs away, what would stand out in 
the solar neighborhood? Undoubtedly, the complex of blue 
star clusters inside the Local Arm, where the Ori OBI and 
Sco OB2 associations, that form the main structure of the 
GB, would just be accessory elements related to the nucleus 
of the complex and its periphery, respectively. 



4 CONCLUSIONS 

From a sample of star clusters younger than 100 Myr, and 
located within a radius of 1 kpc around the Sun, we have 
analyzed the spatial and kinematical structure of the GB. 
The comparative analysis between the populations of stellar 
clusters and OB stars in the GB, indicates that the Ori OBI 
and Sco OB2 associations present a significant difference in 
the number of clusters related to them, and that this differ- 
ence is even more conspicuous in the case of the transient 
stellar condensations. While Ori OBI can be characterized 
as the stellar population associated with the core of a star 
cluster complex (e.g., 30 Doradus), the Sco-Cen complex ap- 
parently shows a star formation mode where the generation 
of isolated stars is dominant. 

In the light of this scenario of hierarchical star forma- 
tion, this difference in the content of stellar clusters must 
have been caused by different physical conditions of the 
primeval clouds and/or inhomogeneities of the gravitational 
potential. The main physical characteristics of these two 
large associations are: 

(i) Different height over the Galactic plane 

(ii) Different content of steUar clusters 

(iii) Different residual velocity vectors 

All these differences can be explained, at least qualitatively, 
by the different position of these two associations with re- 
spect to the main loci of the young stellar clusters that define 
the Orion Arm. According to this scenario, the GB can be 
considered as a partial and biased vision of a much larger 
scale process of stellar formation, which is currently visible 
as a star cluster complex in the region of the Orion Arm 
that is closest to us. 
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Figure 1. Cluster radius vs. age. Tlie filled circles represent our 
probable bound clusters in the age range from 10 to 100 Myr. 
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Figure 2. Pscudo concentration parameter, log(Rc;/Rco) vs. 
cluster radius. Filled circles represent our probable bound clus- 
ters. 
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Figure 3. Distribution of open clusters on the sky (galactic co- 
ordinates). Open and filled circles represent clusters classified as 
LGD and GB members, respectively. Small circles stand for the 
0-B6 GB stars from Paper I. Red crosses and green triangles 
mark, respe ctively, the stars belong ing to Sco OB2 and Ori OBI, 
according to lde Zeeuw et al] l ll999l i. Top panel shows all the clus- 
ters in our sample, whereas bottom panel displays only the prob- 
able bound clusters. 
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Figure 4. Spatial distribution of the sample clusters in the XY 
plane. Symbols as in Figure 3. 
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Figure 5. Spatial distribution of GB clusters (filled circles) in 
the XZ plane. Other symbols as in Figure 3. 
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Figure 6. Same as Figure 4, but including only GB clusters (filled 
circles) with a membership probability higher than 75%. Other 
symbols as in Figure 3. 
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Figure 7. Velocity distribution of GB clusters (filled squares) 
against that of the GB stars (isodensity lines) from Paper I. 
Crosses and open triangles mark, respec tively, the stars belong - 
ing to Sco 0B2 and Ori OBI, according to lde Zeeuw et aI.llll999^ . 
Top panel represents the totality of GB clusters, whereas bottom 
panel displays only those classified as probable transient stellar 
condensations. 
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Figure 8. Distribution of the residual velocities, corrected from 
solar motion and differential rotation, of GB clusters (filled 
squares) against that of the GB stars (isodensity lines) from Pa- 
per I. Crosses and open triangles mark, respec tively, the stars 
belong ing to Sco OB2 and Ori OBI, according to lde Zeeuw et all 
lll999l') . 
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Figure 9. Density map for the star clusters younger than 10 Myr, 
within a box 4 kpc of side centered in the Sun (red dotted circle). 
The Ori OBI and Sco OB2 associations have been superimposed 
on the map, along with their respective residual velocity vectors 
(black lines). 
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Figure 10. Density distribution of the Cluster Formation Index 
(CFI) in the GB plane. The open circles represent the centroids 
of the Ori OBI and Sco 0B2 associations. The solid contour line 
represents the CFI value 0.5, and the dashed contour lines are 
separated by a value of the CFI of 0.05 
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Figure 12. "Reduced GB", showing the system as a schematic 
composition of the centroids of Ori OBI and See OB2 and their 
residual velocities (listed in table [2]| 
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Table 2. Centroids of the OB associations Ori OBI, and Sco OB2, in phase space 



Name 


X 


Y 


Z 


U' 


v 


W 




(pc) 


(pc) 


(pc) 


(kms~^) 




(kms^^) 


Ori OBI 


-369 ± 7 


-178 ± 6 


-132 ±4 


-2±1 


3±1 


1±0.6 


Sco OB2 


118 ±7 


-64 ±6 


29 ±4 


4±1 


-7±1 


0±0.5 



Table 3. Compajrison of the "reduced GB" parameters with their range of values found in the literature 





i 




PQ 






n 


n 


{kms~^kpc~^) 


{kms~^kpc~^) 


Reduced GB 


18 


283 


7 


22 


Literature 


14-27 


271 - 290 


0-29 


12-37 



